Trapping of strangelets in the geomagnetic field 



O 

o 

(N 



L. Paulucc i,i'B J. E. Horvath,2 and G. A. Medina-Tanco^ 

' Instituto de Fisica - Universidade de Sao Paulo 
Rua do Matao, Travessa R, 187, 05508-090, Ctdade Umversitdria 
Sao Paulo SP, Brazil 

^Instituto de Astronomia, Geofisica e Ciencias Atmosfericas - Universidade de Sao Paulo 
Rua do Matao, 1226, 05508-900, Cidade Universitdna 
Sao Paulo SP, Brazil 
^ Instituto de Ciencias Nucleares, Universidad Naeional Autonoma de Mexico 
A. P. 70-543, C.U. Mexico D.F., Mexico 
(Dated: February 5, 2008) 

Strangelets coming from the interstellar medium (ISM) are an interesting target to experiments 
searching for evidence of this hypothetic state of hadronic matter. We entertain the possibility 
of a trapped strangelet population, quite analogous to ordinary nuclei and electron belts. For a 
population of strangelets to be trapped by the geomagnetic field, these incoming particles would 
have to fulfill certain conditions, namely having magnetic rigidities above the geomagnetic cutoff and 
below a certain threshold for adiabatic motion to hold. We show in this work that, for fully ionized 
strangelets, there is a narrow window for stable trapping. An estimate of the stationary population 
is presented and the dominant loss mechanisms discussed. It is shown that the population would 
be substantially enhanced with respect to the ISM fiux (up to two orders of magnitude) due to 
quasi-stable trapping. 
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I. INTRODUCTION 

In a celebrated paper Witten [1] elaborated on the pos- 
sibility 0, H, 0] that systems composed of an deconfined 
Fermi gas of up, down and strange quarks could have a 
lower energy per baryon than iron, thus being absolutely 
stable. This hypothetical state (strange quark matter) 
could be created by weak interactions introducing the 
massive s quark, if the energy "cost" of the mass is com- 
pensated by the availability of a new Fermi sea associated 
to this extra flavor, thus lowering the Fermi energy of the 
u and d quark seas. 

Previous works have shown Q that this stability may 
be realized for a wide range of parameters of strange 
quark matter (SQM) in bulk on the basis of the MIT 
bag model. Calculations also indicate that SQM can be 
absolutely stable within other models, e.g. shell model 
d, 01, or not stable at all More recently, studies 
have indicated that a paired version of SQM, the CFL 
(color-flavor locked) state seems to be even more favor- 
able energetically than the unp aired SQM, widening the 
stability window [9, 1^, [ll|, [l2| ■ 

For the description of finite size lumps of strange mat- 
ter, (termed strangelets) a few terms have to be added to 
the bulk one in the free energy. A surface term suffices 
for A ^ 10^, while other corrections are relevant for the 
lower masses (see [l^ for a recent review). Large lumps 
will have essentially the same number of quarks of bulk 
matter, with a small depletion of the massive strange 
quark resulting in a net positive charge. This is a fea- 
ture also expected for small chunks 0, [l^, which thus 
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resemble heavy nuclei. 

In spite of theoretical controversies, it is generally 
agreed that the ultimate SQM proof must be provided by 
experiments. The experimental searches of strangelets 
started some 20 years ago and have been reviewed re- 
cently in [3, In addition to direct production 
of strangelets in heavy ion collisions [H, [T^, HI, HI, 
[20l . [21I m, m, cosmic rays may contain primaries 
in this state of matter, which could eventually be de- 
tected directly or as a deposition in terrestrial matter 

m HE [13, 113. 

Several cosmic ray events have been tentatively identi- 
fied in the past as primary strangelets (initially the Cen- 
tauro events and the HECRO-81 experiment [29l. [sol. [3l| ) 
for they present features such as their high penetration in 
the atmosphere, low charge-to-mass ratio and exotic sec- 
ondaries [34I. More recently, at least one event recorded 
from the AMS-01 experiment [s^, a mass spectrome- 
ter aboard the shuttle Discovery during a 10-day flight 
in 1998, is being considered as possible detection of a 
strangelets. While it is tempting to identify the primary 
as a strangelet, the inevitable shell effects complicate the 
analysis and preclude any firm conclusion as yet [33 |. It 
is not clear until today to what extent the anomalous 
events can be originated by ordinary primaries or rather 
forcefully require a truly exotic origin. 

Considering the question of existence of strangelets 
among cosmic ray primaries, a few injection (production) 
scenarios have been considered. Witten originally sug- 
gested the merging of compact stars as a likely site [1|. 
In principle, injection spectra and the total mass in the 
galaxy may be calculated knowing the rate of the events 
and the total ejected mass in each of them. These es- 
timates are subject to some caveats, for example, while 
the number of merging systems has been revised upwards 
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[35| . numerical work has shown that a substantial ejec- 
tion of matter is not guaranteed [s^ in a strange star- 
black hole system, and the situation is unclear in the 
case of a fully relativistic SS-SS system, which has only 
been partially addressed [sj because the calculations had 
other goals. On the other hand, strange matter formation 
in type II supernova [ssl] has been preliminarily explored 
and in these events a small fraction of strange matter may 
be ejected. A numerical analysis has shown that the pos- 
sible quark matter component of cosmic rays primaries 
is compatible [s^ with models in which strangelets are 
ejected in either scenario. 

While an uncertain flux from this "contamination" of 
the ISM is expected [i^l, we would like to discuss in 
this paper another likely site to search for strangelets 
of cosmic origin. Much in the same way heavy nuclei 
are present in the earth's magnetosphere bouncing be- 
tween magnetic mirror points, strangelets could also be- 
come trapped in specific regions of the magnetosphere 
and their number density increased respect to the ISM 
flux, provided some conditions for their capture by plane- 
tary magnetic fields are met. This phenomenon is analo- 
gous to the Van Allen belts, and has been first suggested 
in a former study [4l|. A handful of experiments have 
probed the magnetosphere by measuring the fluxes of 
the so-called "anomalous" cosmic ray nuclei, and may 
already place interesting limits to strangelets as well. 
Overall the existence and nature of exotic primaries is 
an important issue. In addition to former and ongoing 
searches, there will be a mass spectrometer placed at 
the International Space Station, the AMS-02 experiment 
p^ . [4^ , with one of its goals to help the identification of 
this exotic component, of crucial importance in testing 
the validity of the Bodmer-Witten-Terazawa conjecture. 
We substantiate below the strangelet belt idea, discuss 
the main features of this population, and advocate for a 
search of this exotic component at definite sites within 
existing uncertainties based on these calculations. 



II. STATES OF IONIZATION AND 
ELECTRONIC RECOMBINATION OF 
STRANGELETS IN THE ISM 

As is well-known, unpaired (also referred as "normal" 
in this work) SQM in bulk contains light u, d and massive 
s quarks in /3-equilibrium. Because of the depletion of the 
more massive s quark, a small fraction of electrons is also 
present to maintain charge neutrality. On the other hand, 
SQM in a paired CFL state is automatically neutral, since 
the equal number of flavors is enforced by symmetry [43 | . 
Actually, a small positive charge is present because of 
the smaller abundance of s quarks near the surface in 
CFL strangelets 13|. Therefore it is natural that CFL 
strangelets will be surrounded by an electronic cloud in 
order to neutralize its total charge, forming an exotic 
atom. The same happens for normal strange matter if the 
strangelet radius is smaller than the electron Compton 



wavelength, a condition satisfied whenever A -C 10^. 

In the following and throughout the whole analysis 
presented here, the strangelet rest mass will be assumed 
to be eo A ~ (930 x A) MeV, with eo the asymptotic 
value of the energy per baryon of strange quark mat- 
ter. We will not consider the fact that the energy per 
baryon number decreases with A in sophisticated model 
calculations, given that the uncertainties found in other 
parameter choices are expected to be much larger than 
the error associated with this approximation. Also the 
strange quark mass is considered to be ms — 1^0 MeV 
and the coupling gap of CFL strange quark matter, 
A = 100 MeV in this exploratory study. With these as- 
sumptions, the net positive charge of strangelets is given 
approximately by Z = OA A (low baryon number regime) 
in the MIT bag model approach for normal strange mat- 
ter and Z = 0.3 A^/^ for the CFL model. 

Strangelets from whatever astrophysical injection 
event would travel through the interstellar medium and 
become ionized by collisions. A simple analysis to evalu- 
ate the degree of ionization of scmi-rclativistic strangelets 
surrounded by electronic clouds due to these interac- 
tions was performed in a Bohr atom approximation. 
Strangelets are partly neutralized by electrons from the 
excitation of the vacuum if Z ^ 100 [i^, but for all 
cases of interest in this work the baryon number range is 
such that we do not have to deal with this effect. 

We considered a two-body collision (incident electron 
- electron in the strangelet cloud) instead of a multibody 
problem, which would be much more difficult to handle. 
The stripping interactions are mainly due to electrons 
with a Maxwellian speed distribution at a temperature 
of ~ 100 K, an average condition of electrons in the ISM. 

The results are shown in figure [1] for strangelets with 
total energy of 1 GeV/A. Considering the average density 
in the interstellar medium to be 1 particle / cw? , the mean 
free path for an electronic collision which may or may 
not result in ionization is of the order of 10^^ cm, which 
is very short on astronomical standards. 

The ionization degree became stable within a travelled 
distance of a few pc for 1 GeV/A strangelets. For ultra- 
rclativistic strangelets (i.e., of the type of candidates that 
would produce a Centauro event [4^ E/A TeV) the 
calculations indicate always full ionization. Furthermore, 
according to the model proposed by Werner and Salpeter 
for the radiation flux in the ISM, the influence of the 
radiation field on ionization of strangelets will be neg- 
ligible unless the strangelet trajectory crosses a region 
containing very energetic photons (i.e. the surroundings 
of a Wolf-Rayet, O and B stars and/or regions of stellar 
formation) . 

We acknowledge that a Bohr atom treatment is a 
crude approach for the electron distribution around the 
strangelet, since it does not include quantum corrections 
as important as the spin-orbit coupling and non-local ef- 
fects, nor relativistic corrections for many electrons bod- 
ies [Z > 40). There is no general expression for these 
corrections applicable in the case of atoms with many 
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FIG. 1: Strangelet effective charge (IGeV/A) versus the 
baryon number A for normal (upper panel) and CFL (lower 
panel) strange matter after interaction with electrons in the 
interstellar medium. 



electrons though; the existent models (e.g., Hartree-Fock 
calculations) are restricted to atoms with few electrons, 
the same happening for experimental corrections. In this 
way, the calculations presented here are rough estimates, 
showing the general trend of the effects rather than pro- 
viding precise numerical values. 

For low-energy particles the electronic capture can be 
as important as the ionization process thus far discussed. 
An approximate cross section for the capture of electrons 
of velocity by a charged particle of atomic number Z 
is given as (ill [H 
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X 6.65 X 10-25 cm^. 



(1) 



where hv k, for Ef, >> I, I and Ee being the electron 
energies while bound to the nucleus and free in the ISM, 
respectively, and irieC^ is the electron rest mass. This 
form of the cross-section for radiative recombination is 
obtained relating the capture of a bare nucleus of charge 
Ze with the capture into the corresponding state of a 
hydrogen atom, which is proportional to the energy of 



the gamma emitted in the process and also to the cross- 
section for the absorption of a quantum of frequency v 
by a ion resulting in emission of an electron of ve- 
locity V, and inversely proportional to the momentum of 
the electron absorbed. In the case of a partially screened 
nucleus, the cross-section is still given approximately by 
equation ([1]), though a special calculation must be per- 
formed to obtain the cross-section for capture into an 
orbital with quantum number uq, usually given in tables 
for ordinary nuclei. 

The "atom" or "ion" formed by capturing an electron 
may also lose this electron in further interactions. For 
light materials the cross section for electron loss can be 
approximately expressed for v > vq [S^I as 



(2) 



where ao = h?' /me^ = 0.53 x 10 ® cm is the Bohr radius 
and Wo = e^/Ti, whereas for intermediate Z materials 
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(3) 



because of the screening effect. 

In summary, a comparison of eqs.(l), (2) and (3) shows 
that electronic capture would only be important for high 
Z strangelets, precisely where this simple picture can no 
longer be applied due to vacuum excitation effects. That 
corresponds to a region in baryon number which we be- 
lieve to be of minimum relevance to the trapped popula- 
tion. 

In summary, these results indicate that we can assume 
total ionization as a good approximation to incoming 
ISM strangelets that could form an ionization belt in the 
magnetosphere. 



III. CAPTURE OF STRANGELETS IN THE 
GEOMAGNETIC FIELD 

The motion of ionized strangelets in the earth magne- 
tosphere can be studied by applying the Stormer theory 
in a dipolar magnetic field. The movement analysis can 
be made in terms of the geomagnetic latitude and the L 
parameter, where L is the equatorial distance of a field 
line to the axis of the dipole measured in units of the 
earth radius. 

The geomagnetic field is not a pure dipole field. In- 
stead, most magnetic models used for studying it include 
nearly 50 terms for describing the potential field from 
which the magnetic field is obtained in a sum of Leg- 
endre functions multiplied by oscillatory coefficients in 
the azimuthal variable. Since the potential field has a 
7.-("+i) dependence, the importance of high-order terms 
decreases rapidly as one moves away from the earth sur- 
face. In this way, the n = 1 term, i. e., the dipole term, is 
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the lowest but dominant term, and most features of the 
trapped radiation theory are analyzed based on a dipole 
field. 

Charged particles with energy of order of MeV in the 
inner part of the magnetosphere {L <^ 10) rotate with a 
much higher frequency than that of typical geomagnetic 
field variation (which varies in time scales of, at most, few 
minutes). Under these conditions, the magnetic moment 
is a conserved quantity (adiabatic invariant). Therefore 
particles with high enough magnetic moments become 
trapped in the dipolar field lines of the geomagnetic field, 
with mirror points placed near the earth poles. 

Particles with mirror points that allow penetration 
in the earth atmosphere can be lost via collisions with 
atoms. All the particles with mirror points placed inside 
the earth radius are obviously lost, meaning that parti- 



cles with laenl < OiE or Ivr 



ae 



< aE, where aeq is the 



equatorial pitch angle, are inside the earth loss cone. 

We will consider collisions mainly with the neutral ni- 
trogen molecule (iV2). The probability of interaction of 
trapped particles which penetrate the atmosphere (suf- 
fering collisions losses) can be taken as 



P{s) - 1 



-s/X{s) 



(4) 



at a certain point s, since each process is probabilistically 
independent, being A the particle mean free path. Gener- 
alizing the previous equation, it is necessary to integrate 
over the particle path. Assuming that all the strangelets 
which collide with particles in the atmosphere are even- 
tually removed from the trapped flux, we express the 
escape probability as 



— J" cr n{s') + s' ds' 



(5) 



where ds = L Re cos sin XdX is the arc along a 

field line, a is the particle cross-section and n{s) is the 
density of particles in the atmosphere at a certain point 
s of the strangelet's path. Since strangelets are hadrons 
we may take their relevant interaction cross-section to be 
geometrical (cr oc A^^^). 

The calculated loss cone for strangelets, assuming an 
exponential profile of the atmospheric density is shown 
in figure [5] for different L reflecting collisions with at- 
mospheric particles and the non-existence of a suitable 
mirror point. It indicates, as expected, that the smaller 
the equatorial pitch angle the easiest it is to remove a 
trapped particle. 

In order for a particle to penetrate a certain region in 
the magnetosphere its energy must be enough to over- 
come the local geomagnetic cutoff rigidity. A solution 
was found by Stormer [5l| describing a special case of 
what he called the "forbidden cone" , which gives the ge- 
omagnetic cutoff rigidity. In this way, the condition a 




FIG. 2: Loss cone for strangelets in the geomagnetic field, i. 
e., particles with mirror point placed at geomagnetic latitudes 
above the square dots (coming from the analysis of collisions 
with atmospheric particles) are eventually removed from the 
trapped population. The full line represent mirror points at 
the earth surface in the dipole field approximation. 



particle must fulfill to have access to a given region of 
the magnetosphere can be written as [s^ 



Rparticie ^ 



59.6 cos* A 



L2[1 + (1--C0S7C0S3 A)i/2]2 



GV (6) 



where A is the latitude and 7 the arrival direction of the 
particle (east - west). 

In the analysis of charged particles trapped in a mag- 
netic field it is usually considered that the motion of a 
given particle is a composition of three different motions: 
the bouncing motion of a guiding center along the mag- 
netic field line; the rotational motion of the particle itself 
around that guiding center; and the longitudinal drift 
of the guiding center. In this way, the condition for a 
triply — adiabatic motion is that the magnetic field in- 
tensity must vary very slowly around a cyclotron orbit, 
imposing a maximum energy allowed for stable trapping. 
The condition that must be imposed for the cyclotron 
radius at the equator is given by 
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Figures[3]and[4]show those bounds for normal and CFL 
strangelets, respectively, for L = 2 in addition to the 
minimum baryon number which is required for strangelet 
stability Q . The existence of a minimum baryon num- 
ber is expected in all models of SQM because the energy 
needed for producing the system increases as the baryon 
number decreases, till it reaches a value above which the 
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strange matter is unstable. The value adopted has been 
■^min = 30 (shown with a vertical line) and may be triv- 
ially altered for any other figure. Strangelets with very 
high baryon number, though allowed for stable trapping, 
arc not likely to be statistically significant for detection 
in the magnctosphere due to a substantial decrease of the 
interstellar flux expected as the baryon number increases. 

The upper bound ([7]) has been enforced in our calcula- 
tions to a 10% confidence level according to observations 
of anomalous cosmic rays L-shell distributions [Esll , and 
we considered E± ^ E, which means we are actually 
underestimating the number of particles that could be 
stably trapped in the geomagnetic field. Obviously, the 
geomagnetic cutoff curve must be below the adiabaticity 
criteria for stable trapping to occur. This is not the case 
for small latitudes, but there is a narrow "window" in 
latitude starting slightly above 30 degrees at L = 2 for 
strangelets from the ISM flux to fulflU the conditions of 
capture and accumulate in regions labelled by the L pa- 
rameter. However, the number of accumulated particles 
is still interesting, as shown in the next section. 

When this calculation is repeated for the case of CFL 
strangelets, the region allowed for stable trapping for 
CFL strangelets has a different shape than that for nor- 
mal strangelets. This feature is due to the strong de- 
pendence of the charge upon A of normal strangelets 
{Z cx A) resulting in constant values if one considers 
momentum per baryon number, whereas for the CFL 
strangelets charge is almost independent of A, leading 
to a ~ A~'^'^ dependence in the momentum per baryon 
number variable. 



A. Trapped strangelet population 

Even though strangelets can be captured and trapped 
in the earth's magnetic field, we must evaluate the pos- 
sible maintenance of a strangelet population to check 
whether there is an increase of the fiux. For this pur- 
pose, we must consider losses mechanisms. 

In addition to the already analyzed losses by collisions 
with neutral atmospheric particles, we have considered 
the inward drift driven by asymmetric fiuctuations of the 
geomagnetic field as a dominant mechanism to diminish 
the strangelet population. 

We will not consider in this work direct pitch angle 
diffusion. Because of their large mass, strangelets are 
less likely to be scattered appreciably in pitch angle by 
collisions. The net result of multiple collisions with atmo- 
spheric particles would be a reduction in the strangelets 
kinetic energy to thermal values and minor changes in 
their pitch angle. Since we have already considered that 
particles bouncing at a radial distance from the surface 
of the earth below the atmosphere height scale (derived 
in section II) would be eventually removed, we are in fact 
replacing a diffusion equation in the cos(aeg) variable for 
a constant loss term (a sink function) directly related, 
though not formally assigned, to pitch angle diffusion. 



Radial diffusion must proceed by fluctuations in the 
third invariant </), which is proportional to L~^, due to 
changes in the electric or magnetic fields that are more 
rapid than the particle drift frequency. Because the gyra- 
tion and bounce periods are much shorter than the drift 
period, the first and second adiabatic invariants are less 
likely to be affected by many of these field perturbations. 

Guided by the existing calculations and observations 
for anomalous cosmic ray nuclei (hereafter ACR) trap- 
ping, we have considered third invariant diffusion due 
asymmetric fiuctuations in the geomagnetic field, which 
is mainly driven by the solar wind pressure (sudden com- 
pression and slow relaxation of the geomagnetic field). 

The diffusion coefficient Dll is determined theoreti- 
cally by taking two consecutive steps [sij. First, one has 
to evaluate the radial displacement suffered by a particle 
under the influence of the field disturbance, which is an 
idealized model of the real disturbances occuring in the 
geomagnetic field. The following procedure is taken in 
order to obtain the diffusion coefficient as a function of 
the statistical features of the disturbances alone. It con- 
sists of squaring this displacement and taking the average 
over several disturbances randomly occurring in time and 
over all possible particle's initial longitudes. 

The diffusion coefficient due to magnetic field fiuctua- 
tions for equatorially trapped particles, with the assump- 
tion of efficient phase mixing [54] can be expressed as 

2 . . 2 p2 r 10 
I5fL = y(y) ^^-L/*i^A(..„/*) (8) 

where Pa{i^) is the power spectral density of the field 
variation evaluated at the drift frequency. For off- 
equatorial particles, the diffusion coefficient presents an 
exponential decay with latitude. 

Already in the case of nuclei, it is known that the com- 
plex geometry and inhomogeneities in the geomagnetic 
field make quantitative calculations ambiguous. The ob- 
served values of the diffusion coefficient and their L de- 
pendence will change with global magnetic activity, and 
magnetic disturbances are known to vary appreciably 
with time. We have assumed a dependence of the 
power spectral density for simplicity [54|. The loss of 
more detailed information associated with this approx- 
imation is that the diffusion coefficient becomes inde- 
pendent of the energy of the particle entering the ge- 
omagnetic field. In this case, the diffusion coefficient 
have a strong dependence on the Mcllwain parameter 
{Dll oc L^°) [bo]. This indicates that its infiuence is 
very important for particles trapped at higher L-shells. 

Typical values for changes in the trapped population 
distribution ranges from a few hours at i = 6 to hundreds 
of days at L — 2. Therefore if strangelets are captured 
by the geomagnetic field their density must be higher for 
lower values of the L parameter, which may result in a 
substantial increase of this population compared to the 
ISM flux. 
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FIG. 3: Restriction curves ([6} and (O for L = 2 in the baryon number vs. momentum plane for normal strangelets incident 
from the east (7 = tt) for different incident directions (reminding that field lines at L = 2 penetrate the earth surface at A = 45° 
in the dipole model). 



Some other losses mechanisms are of less importance 
in short time scales, but have influence on long time 
scales, thus lowering the residence time for trapped par- 
ticles. This includes electrical drift-resonant interactions 
between particles and fields, especially in the pulsation 
frequency or VLF range [55|. Those phenomena are 
highly affected by the solar wind activity. 

The diffusion equation has been employed to study the 
trapped strangelet flux 



df{fi,J,L) _ d 
dt dL 



\Dll d 
dL 



(LV(/i,J,^)) 



(9) 



where / is the distribution function, D^l is given by 
equation ([8]) and /i and J are the adiabatic invariants 
magnetic moment and integral invariant, respectively. 
The relation between the distribution function and the 
flux may be given by j{E,a) = p'^L'^ f{fi, J, L). A sta- 
tionary population requires df /dt — 0, i. e., the assump- 



J 



tion that the source and loss terms are instantaneously 
balanced is valid. 

We assume a steady strangelet injection from the in- 
terstellar medium at L = 6 (the position of the maximum 
distribution function is very insensitive to the chosen L- 
shell parameter for this boundary condition) and derive 
the distribution function shape between this maximum 
and L w 1.05 where it is null (atmosphere particle inter- 
action height), shown in Figure O We are not considering 
diffusion in pitch angle due to interaction of particles with 
electromagnetic waves caused by field variations, which 
alters the first adiabatic invariant. 

The calculations were carried on with two values of the 
flux from the ISM reaching the outer magnetosphere. 

The first one, which will be called "standard", is the 
one that assumes the standard cosmic ray dependence on 
the strangelet flux, E~^-^. The total ISM strangelet flux 
that reaches the earth as estimated by Madsen [40| for a 
binary strange star system coalescence scenario is given 

by 



F w 2 X 10^ m-^ 



yr 



A-"-^'^'Z-^-'max[RsM,RGcr'''A 



(10) 
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FIG. 4: The same as Figure [3] for CFL strangelets. 



where Rsm and Rgc are the solar modulation and 
geomagnetic cutoff rigidities, respectively, and A is an 
uncertain parameter assumed to be of 0(1). In this way, 
the whole flux is fitted with a E~^-^ dependence with the 
constraints of minimum and maximum energy respecting 
the values i?™,„ = bMVA/Z and Rmax = IQ^GV 

The second calculation, which will be called "im- 
proved" , considers a more detailed characterization of the 
differential flux, where for the region of interest in this 
work (rigidities of few GV), the strangelet flux actually 
increases with a slope of i?^'^. This flux was obtained 
from a flt to reference ,40]. 

In either way, the flux entering the region of the mag- 
netosphere at L„iax has to fulfill the restrictions imposed 
for stable trapping in the pitch angle and geomagnetic 
latitude of incidence 



Fz„= / d\P{X) / daeqP{a,g)xF (11) 



The efficiency factors, P{\) and P{aeq) may be easily 
identified: P(A) gives the fractional area of the spherical 
section suitable for trapping discussed previously 



^ 2L^{-cos\)f^-dcp 

where the factor 2 comes from the symmetry in 9 for 
both hemispheres (north, south). P{aeq) limits the num- 
ber of particles entering the specific region of the mag- 
netosphere with an appropriate pitch angle to avoid the 
loss cone as already discussed. We have also assumed an 
isotropic fiux of particles, since there is no theoretical pre- 
diction pointing to any anisotropy in the arrival direction 
of strangelets, which means that jo{cosaeq) = constant 
is a reasonable hypothesis: 

P{a,q) = ^—j^ 

Jo ^eg ^OL^q 

where the factor 4 stands for the symmetry in the con- 
dition for a given particle to belong to the loss cone: 

Q^eq I ^losscone and [tT tteg | ^loss cone- 

Solving the differential equation ^ and obtaining the 
corresponding flux inward (in the -e^. direction) for every 
L, it is possible to determine the mean particle density 
at a given shell and, in this way, the trapped strangelet 
flux. 
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FIG. 5; The distribution function (in arbitrary units) for 
strangelets trapped in the geomagnetic field as a function of 
L is obtained as the solution for the difi^erential equation |9] 
with the boundary conditions f{Lmax) given by the incom- 
ing flux (see text for details) and f{Lmin) = (correspond- 
ing to the height scale of interaction with atmospheric par- 
ticles). The position of the peak (around L = 1.3) does not 
change with the change in the A/Z relation (CFL and non- 
CFL strangelets) nor with a change in the energy and A of 
the strangelets. 





L = 1.3 


L = 2 


Normal 


1.28 X 10"^^ 


4.34 X 10"^'' 


CFL 


3.95 X 10"^" 


1.34 X 10~^^ 



TABLE I: Mean particle flux in units of 
partcm~^ s"^ sr~^ {MeV/A)~^ for strangelet stationary 
population at i = 1.3 and L = 2 calculated with the 
"standard" flux. 



The results are summarized in Tables 1 and 2 for 
L = 2 (ACR belt location) and L ~ 1.3 (location of 
the maximum of the distribution function) for the exam- 
ple of strangelets of ^ = 100 and energy corresponding 
toR= IGV. 

The position of the peak of the distribution function in 
the geomagnetic field (around L = 1.3) is quite robust, it 





L = 1.3 


L = 2 


Normal 


3.83 X 10-^'' 


1.3 X 10"^'^ 


CFL 


1.64 X 10"" 


5.55 X lO"^'^ 



TABLE II: Mean particle flux in units of 
partcm~^ sr~^ {MeV/A)~^ for strangelet stationary 
population at L = 1.3 and L = 2 calculated with the 
"improved" flux. 



does not appreciably change with the change in the A/Z 
relation (CFL and normal strangelets), nor with a change 
in the energy and baryon number of the strangelets. This 
could be a consequence of the assumption of the power 
spectral density as being proportional to v^^, what ren- 
ders the diffusion coefficient independent of the particle 
energy, therefore modifying the energy of the particles 
does not affect their diffusion properties. 

We observe that the trapped population is slightly 
more favored if strange quark matter is in the CFL state, 
the difference between the trapped fluxes for the two 
species increases with decreasing energy exponent in the 
incident flux. It happens due to the dependence on the 
baryonic number of the interstellar flux of strangelets (eq. 
[TO]). Since the rigidity interval for stable trapping is the 
same for both states for it only depends on the geometri- 
cal characteristics of the geomagnetic field, the difference 
on the number of particles trapped strongly depends on 
the difference in the incoming flux. This dependence of 
the integrated flux on the number of baryons that can 
be expressed as Fjsm oc (0.125)~^-^ ^-i.667 Fjsm « 
(0.3)"^'^ ^-1.267 £pj. normal and CFL strangelets, respec- 
tively. In this way the flux of paired CFL strangelets is 
lower than those without pairing, but only for low baryon 
number {A <~ 13), that is, in a region where it is be- 
lieved strangelets are not stable at all. In the stability 
region the flux for CFL strangelets is always higher than 
for normal strangelets resulting in a higher trapped den- 
sity. In this way, the smaller difference seeing between 
strangelets with and without pairing for the improved 
flux when comparing to that for the standard flux is ex- 
plained by the smaller difference in the incoming flux due 
to the softer dependence on the atomic number (the stan- 
dard flux depends on E~'^-^, which for the same rigidity 
depends on the particle's atomic number; instead, the 
improved flux depends on R^'^ and the analysis was per- 
formed in terms of same rigidity). 

Additional considerations are relevant for the fate of a 
trapped population of strangelets. It is well-known that 
the solar wind has a strong influence on the ACR flux 
upon the earth. The most abundant ACR heavy ion, 
oxygen, shows a strong intensity variation with the solar 
cycle, having its interstellar flux of 8 — 27 MeV/nucleon 
lowered up to two orders of magnitude during periods of 
solar maximum activity [s^ . During solar minimum, the 
trapped flux at the earth magnetosphere is of the order 
of ^ 5 X 10^^ particles cm~^ sr~^ (MeV/nucleon)~\ 
corresponding to an enhancement factor of ^ 15 [57| . 
this experimental value being somewhat below the the- 
oretical expected one (higher than 25 [1^). The oxy- 
gen component corresponds to about 80% of the trapped 
ACR, while the C/0, N/0 and Ne/O abundance ratio 
are < 0.005, - 0.10-0.15 and - 0.02-0.03, respectively. 

With the results obtained in this study, the trapped 
flux of strangelets at L < 2 would be of or- 
der 10^^'^ -IQ-^^ particles 0771^"^ sr~^ s~^ {MeV/Ay'^ at 
rigidity R = 1 GV for strangelets of baryon number 
A = 100. This represents an enhancement factor for 



9 



trapped flux in the regime of steady-state population 
comparing to the interstellar flux at the same energy and 
A of order 10 and 10^ for strangelets trapped at L = 2 
and L — 1.3, respectively, the values for CFL strangelets 
being of about twice the one for normal strangelets 
{q ~ 5.5 and 11, and q ^ 162 and 314 for CFL and nor- 
mal strangelets at L = 2 and L = 1.3, respectively). This 
results show that the strangelet flux could be as high as 
a factor 10000 lower than that expected for carbon dur- 
ing periods of maximum solar activity. Although we did 
not consider the solar modulation in our calculations, it 
would act significantly over those low energy strangelets 
[4^, the region of interest in this study. In this man- 
ner, it could have an important influence, similar to that 
detected for oxygen, on the trapped density. 

The advantage of a search for trapped strangelets in 
the geomagnetic fleld performed during the solar maxi- 
mum activity whether they are an important component 
of the radiation belt or are to be measured penetrating 
the atmosphere towards to surface of the earth would be 
the reduced component of ACR, which could reduce dead 
time losses in the detectors and possibly render a clearer 
identification of the primaries. 

The proposed and widely accepted model for ACR 
trapping [53| assumes that the high mass-to-charge ra- 
tio of singly-ionized ACRs enables them to penetrate 
deeply into the magnetosphere. ACRs with trajectories 
near a low altitude mirror point interact with particles 
in the upper atmosphere, loosing one or all their remain- 
ing electrons. After stripping, the particle gyroradius is 
reduced by a factor of 1/Z, and the ion can become sta- 
bly trapped. As stated above, the results presented here 
were obtained assuming fully- ionized strangelets, which 
have just the "right" features to become trapped. How- 
ever, some fraction of the strangelets should reach earth's 
atmosphere with an effective charge slightly below their 
atomic number and suffer a process of interaction simi- 
lar to ACR's, which is much less dependent on the pitch 
angle and other variables. Finally, there is also the possi- 
bility of quasi-stable trapping of ions with energies high 
enough not to obey condition ([7]), but not too high as to 
penetrate the magnetosphere without suffering any sig- 
nificant depletion in their incident direction. These two 
additional mechanisms could result in a further increase 
in the number of trapped strangelets. 

Conclusions 

From the analysis presented here we conclude that non- 
relativistic strangelets with A <^ 10^ already ionized 
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by collisions with electrons in the ISM could be stably 
trapped by the geomagnetic field. Assuming the exis- 
tence of a strangelet contamination in the ISM, its in- 
jection in the solar system and given the geomagnetic 
geometry and the interaction of the magnetic field with 
the solar wind, it looks very likely to have this radiation 
belt surrounding the planet. If strangelets are to be a 
component of the anomalous cosmic ray belt at i ^ 2, 
we have shown that, even considering the approximations 
taken during the calculations presented here (which have 
the main consequence of averaging the trapped popula- 
tion's behavior), those particles would be present with 
an enhancement factor comparing with the interstellar 
fiux of order 10^ and if we consider a new particle belt 
(a strangelet belt) at i ^ 1.3, the enhancement factor 
could be as high as order 10^ in a stationary popula- 
tion scenario |6l| . These exotic baryons could in princi- 
ple be detectable in the earth magnetosphere depend- 
ing on the chosen parameters for each of the experi- 
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Overall, we believe our estimates to be on the con- 
servative side of the trapped flux, making the search of 
trapped strangelets a feasible but difficult task. Needless 
to say, the detection of those trapped particles having low 
Z/A ratio would be extremely important for determining 
the properties of cold, dense baryonic matter. 



Acknowledgments 

L.P. acknowledges the financial support received from 
the Fundagao de Amparo a Pesquisa do Estado de Sao 
Paulo. J.E.H. wishes to acknowledge the CNPq Agency 
(Brazil) for partial financial support. 



[4] H. Terazawa, Tokyo U. Report, pp. INS-336 (1979). 
[5] E. Farhi and L. Jaffe, Phys. Rev. D 30, 2379 (1984). 
[6] E. P. Gilson and R. L. Jaffe, Phys. Rev. Lett. 71, 332 
(1993). 



10 



[7] J. Madsen, astro-ph/9809032. 
[8] M. Buballa, hep-ph/0402234. 

[9] G. Lugones and J. Horvath, Phys. Rev. D 66, 074017 
(2002). 

[10] M. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B 

537, 433 (1999). 
[11] R. Rapp et al., Ann. Phys. (N. Y.) 280, 35 (2000). 
[12] K. Rajagopal and F. Wilczeck, for an overview of the 

CFL state, hep-ph/0011333. 
[13] J. Madsen, J.Phys. G 28, 1737 (2002). 
[14] R. Khngenberg, J. Phys. G 27, 475 (2001). 
[15] E. Finch, nucl-ex/0605010. 

[16] J. Thomas and P. Jacobs, A guide to the high energy ion 

experiments, UCRL-ID-119181. 
[17] A. Rusek et aL (for the E886 collaboration), Phys. Rev. 

C 54, R15 (1996). 
[18] G. Van Buren (for the E864 collaboration), J. Phys. G 

25, 411 (1999). 

[19] J. Belz et al. (for the BNL E888 collaboration), Phys. 
Rev. Lett. 76, 3277 (1996). 

[20] F. Dittus et al. (for the NA52 collaboration), in Inter- 
national Conference on Strangeness in Hadronic Matter, 
edited by J. Rafelski (American Institute of Physics, New 
York (AIP 340), 1995), p. 24. 

[21] G. Appelquist et al., Phys. Rev. Lett. 76, 3907 (1996). 

[22] G. Ambrosini et al., Nucl. Phys. A610, 306c (1996). 

[23] R. Khngenberg, J. Phys. G 25, R273 (1999). 

[24] R. Arsenescu et al.. New J. Phys. 4, 96 (2002). 

[25] A. De Rujula and S. L. Glashow, Nature 312, 734 (1984). 

[26] M. Bruegger et al.. Nature 337, 434 (1989). 

[27] M. C. P. Isaac et al., astro-ph/9806147. 

[28] Z.-T. Lu et al., nucl-ex/0402015. 

[29] J. D. Bjorken and L. McLerran, Phys. Rev. D 20, 2353 
(1979). 

[30] M. Rybczynski, Z. Wlodarczyk, and G. Wilk, hep- 
ph/0109225. 

[31] T. Saito et al., Phys. Rev. Lett. 65, 2094 (1990). 
[32] M. Rybczynski, Z. Wlodarczyk, and G. Wilk, hep- 
ph/0410065. 

[33] V. Choutko (for the AMS-01 Collaboration), in Proc. 
28th Internat. Cosmic Ray Conf., Tsukuba, Japan (Uni- 
versal Academic Press, Tokyo, 2003), p. 1765. 

[34] J. E. Horvath and L. Paulucci, J. Phys. G 32, B13 (2006). 

[35] V. Kalogera et al., Astrophys. J. Lett. 601, L179 (2004). 

[36] W. Kluzniak and W. Lee, MNRAS 335, L29 (2002). 

[37] F. Limousin, D. Gondek-Rosinska, and E. Gourgoulhon, 
gr-qc/0411127. 

[38] O. G. Benvenuto and J. E. Horvath, Phys. Rev. Lett. 63, 
716 (1989). 



[39] G. A. Medina- Tanco and J. E. Horvath, Astrophys. J. 

464, 364 (1996). 
[40] J. Madsen, Phys. Rev. D 71, 014026 (2005). 
[41] H. C. Rosu, hep-ph/9410028. 

[42] Home page of AMS experiment, URL 

http : / /ams . cern . ch 
[43] J. Madsen, hep-ph/0111417. 

[44] K. Rajagopal and F. Wilczek, Phys. Rev. Lett. 86, 3492 
(2001). 

[45] J. Madsen and J. M. Larsen, Phys. Rev. Lett. 90, 121102 
(2003). 

[46] C. Lattes and K. Fujimoto, Phys. Repts. 65, 151 (1980). 
[47] M. W. Werner and E. E. Salpeter, MNRAS 145, 249 
(1969). 

[48] H. S. Massey and E. H. S. Burhop, Electronic and lomc 
Impact Phenomena (Oxford University Press, 1952). 

[49] H. A. Bethe and E. E. Salpeter, Handbuch der Physik, 
vol. 35 ( Springer- Verlag Berlin, 1957). 

[50] N. Bohr, Kgl. Dansk. Vid. Selsk. 18, 8 (1948). 

[51] C. Stoermer, The polar aurora (Claredon Press, Oxford, 
1955). 

[52] J. H. Adams, R. Silberberg, and C. H. Tsao, Tech. Rep. 
NRL Memorandum Report 4506, Naval Research Labo- 
ratory, Washington DC 20375-500, USA (1981), See also, 
for example, D. F. Smart and M. A. Shea, Adv. Space 
Res. 14(10), 787 (1994). 

[53] A. J. Tylka, in Proc. 23rd Internat. Cosmic Ray Conf., 
Calgary, Canada (1993), vol. 3, p. 436. 

[54] M. Walt, Introduction to geomagnetically trapped radia- 
tion (Cambridge Press, 1994). 

[55] W. Baumjohann, private communication. 

[56] S. Biswas, Space Science Rev. 75, 423 (1996). 

[57] V. V. Bobrovskaya et al.. Cosmic Research 39, 98 (2001). 

[58] J. B. Blake and L. M. Friesen, in Proc. 15th Internat. 
Cosmic Ray Conf., Plovdiv, Bulgaria (1977), vol. 2, p. 
341. 

[59] C.-G. Falthammar, in Earth's Particles and Fields 
(Ed.B.M. McCormac, Reinhold Book Co., 1968), pp. 
157-169. 

[60] Radial diffusion caused by random variations in the po- 
tential electric fields have a softer dependence on the 
Mcllwain parameter, being the resultant diffusion coeffi- 
cient Dll oc [H. 

[61] It must be noted, however, that for the ACR belt the 
theoretical values for the enhancement factor are some- 
what higher than the experimental ones, a feature that 
would probably hold for the strangelet belt as well. 



